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P ! Abstract 

& ■ 

We report the first result of the clustering analysis of Sloan Digital Sky Survey (SDSS) quasars. We 
compute the two-point correlation function (2PCF) of SDSS quasars in redshift space at 8/i -1 Mpc < s < 
500/i _1 Mpc, with particular attention to its baryonic signature. Our sample consists of 19986 quasars 
extracted from the SDSS Data Release 4 (DR4). The redshift range of the sample is 0.72 < z < 2.24 (the 
mean redshift is z — 1.46) and the reddening-corrected i-band apparent magnitude range is 15.0 < m, irc < 
19.1. Due to the relatively low number density of the quasar sample, the bump in the power spectrum 
due to the baryon density, fib, is not clearly visible. The effect of the baryon density is, however, to 
distort the overall shape of the 2PCF. The degree of distortion makes it an interesting alternate measure 
of the baryonic signature. Assuming a scale- independent linear bias and the spatially flat universe, i.e., 
fib + fid + fiA = 1 , where fid and fiA denote the density parameters of dark matter and the cosmological 
constant, we combine the observed quasar 2PCF and the predicted matter 2PCF to put constraints on fib 
and fiA. 

Our result is fitted as 0.80 — 2.8fib < fiA < 0.90 — 1.4fib at the 2<r confidence level, which is consistent 
with results from other cosmological observations such as WMAP. The "mean" bias parameter of our quasar 
sample is 1.59(78 _1 (for fib = 0.04 and fiA = 0.7), where er 8 is the top-hat mass fluctuation amplitude at 
8/i _1 Mpc. We also estimate the corresponding bias parameter of quasars at z = 0, &Qso,Fry(0), assuming 
Fry's bias evolution model. For fi b = 0.04, fi A = 0.73 and fi d = 0.23, we find fe QS o,Fry(0) = 0.54 + 0.83CT8" 1 
which is valid for 0.6 < erg < 1.0. 

Key words: galaxies: quasars: general — cosmology: large-scale structure of universe — cosmology: 
observations — methods: statistical 
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1. Introduction 

Quasars are intrinsically bright and preferentially located at relatively high rcdshift. They arc, therefore, naturally 
suited for exploring the structure of the distant universe out to redshifts z ~ 3, a regime not easily accessible from 
existing wide-field galaxy surveys. While surveys of the universe up to z ~ 6 have significantly advanced in recent 
years using Lyman break-galaxies and Lyman-a emitters (e.g., Steidel et al. 1999; Shimasaku et al. 2003; Ouchi et al. 
2004a; Ouchi et al. 2004b), they remain limited in terms of the areas they cover (barely extending to O(100)/i _1 Mpc). 
As such they are not well suited to measuring the large scale clustering at high redshift (i.e. they are sample variance 
limited). Quasar samples, in contrast, can be surveyed over very large volumes and, therefore, play an important and 
complementary role in measuring the evolution of structure in the universe on the largest accessible scales. 

The primary goal of this quasar clustering analysis is to directly probe and quantify the large-scale structure of the 
universe on scales up to a few /i _1 Gpc, where h denotes the Hubble constant in units of 100kms _1 Mpc _1 . Eventually 
such analyses may enable us to answer fundamental questions including the homogeneity of the universe on Gpc scales 
or the validity of the cosmological principle (e.g., Lahav, Suto 2004). Also the clustering signal provides important 
clues to the formation and evolution of quasars themselves and may shed light on the growth history of their central 
massive black holes. 

Although there are several seminal papers addressing quasar clustering (Osmcr 1981; Shanks et al. 1987; Boyle, 
Mo 1993; Croom, Shanks 1996), those previous quasar samples were fairly small which limits their statistical power in 
constraining large scale structure. The 2dF QSO Rcdshift Survey (2QZ:Croom et al. 2001; Croom et al. 2002; Outram 
ct al. 2003) and the SDSS are among the largest statistical samples which should significantly advance the statistics 
of quasar clustering analysis. 

This is the first of a scries of SDSS quasar clustering papers and is aimed at exploring the impact of the baryon 
density and the cosmological constant on the quasar 2PCF on scales of ~ 100/i _1 Mpc. In this sense, our current 
approach is quite complementary to the existing galaxy clustering studies by the SDSS collaboration (e.g., Tegmark 
ct al. 2004; Pope et al. 2004). The samples and techniques are complementary: the galaxies and quasars are selected 
using different algorithms, they explore vastly different volumes, the number densities of the objects are very different, 
and they cover very different redshift ranges. It would be interesting to see if both data sets yield consistent results 
on constraints on the cosmological parameters. 

We will address the dependence of the quasar 2PCFs on redshift and luminosity and comparison with those of 
galaxies and AGNs in a separate paper (Connolly et al. 2005). To be specific, we compute the 2PCF of quasars 
from SDSS Data Release 4 (DR4) on scales from 8/i _1 Mpc to 2/i _1 Gpc and perform statistical analyses that pay 
particular attention to the embedded baryonic signature. From recent progress in the theoretical and observational 
understanding of the initial conditions present in the early universe, the 2PCF of the matter is well described. Two 
characteristic scales in the matter 2PCF, £ m (r,z), are present at approximately 100/i _1 Mpc; the baryon bump and the 
zero-crossing point (denoted by s zer o in what follows). These features are sensitive to the baryon density parameter, 
fib, and the matter density parameter, O m , where we define O m = fib + fid, with fid being the density parameter of 
cold dark matter (CDM). 

We derive constraints on fib and fl\ assuming a spatially flat universe (i.e., fib + fid + fiA = 1) throughout this 
paper. Our attempt may be partially interpreted as yet another cosmological test to determine the values of fib and 
fiA independent of existing methods. In this respect, our current results discussed below are not stringent enough to 
improve the previous constraints. Nevertheless we believe that the attempt to look for the baryonic signature in the 
cosmological distribution of quasars is of scientific interest on its own. 

2. SDSS quasar sample 

2.1. The Sloan Digital Sky Survey 

The Sloan Digital Sky Survey (SDSS;York et al. 2000 for a technical summary) provides the largest homogeneous 
sample of quasars currently available, and thus is best suited for our current purpose. The wide-field five band 
photometric system (Gunn ct al. 1998; Fukugita et al. 1996; Hogg et al. 2001; Ivezic et al. 2004; Smith et al. 2002) 
enables the identification of quasars with high efficiency and completeness. The precision astrometry (Pier et al. 2003) 
maintains the accuracy of the clustering analysis even in 100/i _1 kpc scale. The SDSS is also excellent in terms of 
the homogeneity of the spectroscopic sample because its spectroscopic targets are selected consistently on the basis of 
their photometric data. The details of quasar target selection algorithm are described by Richards et al. (2002), and 
the spectroscopic tiling algorithm is described by Blanton ct al. (2003). 

2.2. Quasar sample selection 

The spectroscopic sample of quasars in the DR4 is selected according to several different spectroscopic target selection 
criteria (see Table 27 of Stoughton et al. 2002). We consider only the point sources, and extract our quasar sample 
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for the current analysis on the basis of four conditions; first, to ensure a homogeneity of the sample, we choose those 
spectroscopic quasars whose "primTarget" flag (see Table 6 of Stoughton et al. 2002) indicates either "QSO.CAP" 
or "QSO_SKiRT" (sec Fig.l and section 3.4 of Richards et al. 2002). Second, their i-band magnitude corrected for 
the Galactic reddening should satisfy 15.0 < mi. rc < 19.1, which are the limits used by the quasar target selection 
algorithm. We adopt the reddening correction of Schlegel, Finkbeincr, Davis (1998). Although this correction is 
calculated adopting the typical spectrum energy distribution of elliptical galaxies and the old SDSS filter response 
functions the resulting difference is negligibly small. Third, the redshift limits of this sample are 0.72 < z < 2.24. In 
the SDSS photometric system, the color locus of quasars in the redshift range 2.4 < z < 3.0 crosses the stellar locus and 
thus the survey completeness becomes significantly degraded in that redshift range. This defines the upper redshift 
limit of our sample. The lower limit is chosen so that the absolute magnitude of all the selected quasars should be 
brighter than -23 in i-band. Finally, we exclude the southern sky region because the selection of the region is different 
from that of the northern sky. The total number of selected quasars is 19986. 

The SDSS photometric data have two different versions; "BEST" version denotes the highest quality data at the 
time of the data release, while "TARGET" version is the data at the time the target selection algorithm was run for 
that part of the sky (i.e., the photometric data may be updated at the time of the data release). Our current sample 
uses the "primTarget" flag which is based on the "TARGET" version. In reality, however, the two versions may lead 
to some systematic effect due to the mismatch of the time-dependent selection function. We compared BEST and 
TARGET using the SDSS Data Release 3 (DR3: see Abazajian et al. 2004) QSO catalog (Schneider et al. 2005), 
and found that the difference does not affect the resulting 2PCFs within the current errorbars (see Fig. 4 below). We 
note that our sample in the DR3 region slightly differs from DR3 QSO catalog (Schneider et al. 2005) because we 
have not manually inspected all objects classified by the pipeline as quasars. The number of quasars in a subset of 
our DR4 sample that is located in the DR3 region is 16656, out of which 32 objects are not listed in the catalog of 
Schneider et al. (2005) and 9 objects have incorrect redshifts. This difference of the quasar samples, however, has a 
negligible effect on the estimated 2PCF as will be shown below. 

The equatorial coordinate distribution of our quasar sample is plotted in Figure 1, and a projected wedge diagram 
of the quasars is shown in Figure 2. Figure 3 shows the scatter plot of redshift and i-band absolute magnitude ( Upper 
panel) and the redshift histogram (Lower panel). 

Giving the errorbars in Fig. 4, the photometric completeness, fiber collisions and the effect of errors in the reddening 
correction are not supposed to change our result. Those issues will be fully addressed by Connolly et al. (2005) who 
found that none appears to significantly affect the quasar 2PCFs presented below. 

2.3. Analysis 

We use the Landy-Szalay estimator (Landy, Szalay 1993) in computing the 2PCF from the quasar sample: 
, t . \ QQ{si < s < s i+ i) QR{s l < s < s i+1 ) 

4QSO,obs < s < s l+ i ) = — — — — — 7-2-^7 — 7 + 1, (1) 

RR(s l < s < s l+1 ) RR(Si < s < s l+1 ) 

where QQ, QR, and RR are the number of corresponding pairs of quasars and random particles that arc normalized 
by Nq(Nq — l)/2, NqN r , Nr(Nr — 1)/2, respectively (Nq and N R are the numbers of quasars and random particles). 
The accuracy and convergent properties of the Landy - Szalay estimator are discussed in the appendix of Kayo et al. 
(2004) and in Kerscher, Szapudi & Szalay (2000). 

The redshift distribution of random particles is determined so as to reproduce that of the quasars. The latter is 
estimated by averaging over a finite width of redshift, Az = 0.08. We find that other choices of Az = 0.02, Az = 0.04 and 
Az = 0.10 result in almost identical ^Qso,obs(s). The number of random particles, Nr, is about 600,000, which is found 
to provide robust estimates of ^Qso,obs(s) on the scales s > 8/i _1 Mpc. We evaluate the covariancc matrix, CV, , from 40 
jack-knife re-samplings (see e.g., Lupton 1993). In order to ensure the statistical significance, the comparison between 
the observed and the predicted 2PCFs needs to be done after the finite binning of the pair-separation. Empirically, we 
compute the 2PCFs at 6 different binning offsets; the i-th boundary of the n-th binning offset (n = 1 . . . 6) is defined as 

Si (n) = 10 a3i+0 - 05( ™- 1) /i- 1 Mpc. (2) 

We then compute the averaged 2PCFs of theoretical predictions (described in section 3) as follows: 

J H+1 ^ Q so(s')s ,2 ds' 

Sqso (si < s < s l+1 ) = Sl j-« +Vw • ( 3 ) 

In practice, for a given n-th offset of the binning, we use bins whose boundaries satisfy 7.9ft.~ 1 Mpc < s, < 502/i~ 1 Mpc 
for the following analysis. Figure 4 plots the 2PCFs of our samples for the n = 1 binning offset adopting Sl m = 0.3, 
Q,\ = 0.7 and h = 0.7. We also plot the 2PCFs of a subset of our DR4 sample that is located in the DR3 region, and 
two 2PCFs of BEST and TARGET samples of the DR3 quasar catalog (Schneider et al. 2005) for comparison. The 
four sets of 2PCFs agree with one another to within la. Throughout the paper, we show the 2PCFs in log-log plots. 
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In order to show the negative part of 2PCFs, we plot the range of —10 2 < £ < — 10 5 in the lower panel of each figure 
by taking the logarithm of its absolute value. 

3. Theoretical predictions 

3.1. Parameter dependence of the matter 2PCF 

Since the baryon bump is the counterpart of the acoustic oscillation peaks in a matter power spectrum, the analysis 
using the 2PCF should, in principle, be equivalent to that with the power spectrum. Nevertheless we perform the 
analysis with the 2PCF for the following two reasons. The first advantage is based on the different behavior of the 
baryonic signature in the power spectrum and 2PCF. While the acoustic oscillation peaks in a matter power spectrum 
show up in several overtones, they translate coherently to the baryon bump in the 2PCF roughly at one specific scale. 
Thus the resulting signature is expected to be more prominent in the 2PCF. The second is related to the inevitably 
complicated geometry of the survey volume boundary. Strictly speaking, the estimate of the power spectrum from 
the observed sample has to assume a periodic boundary condition unless resorting to a very sophisticated estimation 
methodology (e.g., Pope et al. 2004). On the other hand, the 2PCF can be reliably estimated by taking account of the 
arbitrary boundary shape of the survey volume in a relatively straightforward fashion. For this purpose, we simply 
distribute random particles over the survey volume according to the selection function of quasars (the lower-panel of 
Fig. 3) and perform standard pair-counting. This is why we adopt the 2PCF in the present analysis. 

In locating the baryon bump and the zero-crossing point, it is important to have accurate predictions of the matter 
2PCF. Here we adopt an empirical fitting formula for the CDM transfer function incorporating the effect of Ob 
(Eisenstein, Hu 1998). The theoretical model of 2PCFs depends on O m and Ob and not on Oa, while the estimate of 
the observed 2PCFs requires the values of O m and Oa. To simplify the analysis, therefore, we adopt a spatially-flat 
universe, O m + Oa = 1, and choose Ob and Oa as two model parameters. 

Figure 5 shows £ m (r,2 = 0) for different values of Ob (Oa is fixed to 0.73). Baryon bumps are always visible around 
100/i -1 Mpc (except for the case of Ob = 0) and their height and position are sensitive to the value of Ob. Figure 6 
also plots £ m (r, z = 0) but for different values of Oa assuming Ob = 0.04. It is clear that s zero is very sensitive to Oa 
and, to a lesser extent, to Ob- 

Figure 7 shows the dependence of £ m (r, z = 0) on h. Interestingly, s zcro is quite insensitive to the value of the Hubble 
constant. If Ob = 0, the CDM linear matter spectrum is approximately given as a universal function of the scaled 
wavenumber, q = fc/(O m /i 2 ) implying the scaling of s zoro [/i _1 Mpc] oc (O m /i) _1 (see also Matsubara 2004 for related 
discussion on the parameter dependence). We made sure that this scaling holds for Ob < 0.01 but is replaced by a 
different empirical relation of s zero [/i _1 Mpc] oc O" 1 = (1 — Oa) -1 for Ob > 0.02 (the second equality comes from our 
assumption of the spatially flat universe) . We note here that the above feature may be very useful in breaking the 
degeneracy in the constraints on O m and h. 

3.2. Predicting 2PCF of quasars on the light-cone 

In order to predict the observable 2PCF of quasars £qso(s) from the matter 2PCF £ m (r,z), one has to take into 
account three major effects, (i) the light-cone effect, (ii) redshift-space distortion and (iii) quasar biasing. Since 
the formulation and the prescription of incorporating these effects have been extensively discussed in the literature 
(Yamamoto, Suto 1999; Suto et al. 1999; Hamana, et al. 2001; Hamana et al. 2001; Lahav, Suto 2004), we simply 
summarize the results below. 

3.2.1. Light- cone effect 

Since the current quasar sample extends from z = 0.72 to z = 2.24, the clustering evolution within the survey volume 
cannot be neglected; for instance, the linear growth rate of density fluctuations, D(z), in the universe with O m = 0.27 
and O a = 0.73 varies from D(0)/Z)(0.72) = 1.41 to D(0)/D(2.24) = 2.49. This implies that the corresponding matter 
2PCF in the linear regime evolves by a factor of [D(0.72)/D(2.24)] 2 w 3.1 between the two edges of the survey volume. 
Although wc could reduce the effect in principle by dividing the entire sample into several subsamplcs in a narrow 
redshift range, it would appreciably degrade the statistical significance of the result. This is why one has to carry out 
the proper averaging procedure in predicting £qso(s) theoretically. To be specific, we follow the formula described by 
Hamana, et al. (2001): 



where z m i n and z max denote the redshift range of the survey, dV c /dz is the differential comoving volume clement, <p(z) 
is the selection function of quasars, no(z) is the number density of quasars, and Cqso(sj^) is the 2PCF of quasars in 




(4) 
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redshift space at z. 

In practice, we rewrite equation (4) as 

dz ' ' I dN QSO[z) 



dz J \ dz 

In the above expression, dNQgo(z)/dz is the redshift number distribution of quasars which is directly estimated from 
the data (Lower panel of Fig. 3), and dV c /dz can be easily computed once a set of cosmological parameters is specified. 

3.2.2. Redshift- space distortion 

The next task is to convert the theoretical predictions in real space to those in redshift space. Since the current 
paper focuses on scales beyond ~ lO/i^Mpc, the nonlinear effect (i.e., fingers of god) is not important and we may 
just consider the Kaiser effect (Kaiser 1987): 

2 flM| 1 2 ^ f _ ^ Qf ^_ 1 dlnD(z) 



^QSo(s,z)={l + -P(z) + -P'(z)j^ QSO (r,z), ^ z) = b QSO (z) d\na ' (6) 

where r and s denote the quasar pair-separations (comoving) in real and redshift spaces, respectively, a is the scale 
factor, and 6qso( z ) is the biasing factor of quasars at z as described below. This expression assumes the distant- 
observer approximation, and can be justified for pair-separations less than 10° on the sky (Matsubara 2000). For 
reference, at the median redshift of z ~ 1.4, the comoving scale subtended by an angle of 10° is w 511/i _1 Mpc in the 
f2 m = 0.3 and Oa = 0.7 model. This is why we only use the bins less than 502 /i _1 Mpc where the distant-observer 
approximation can be reasonably justified. 

3.2.3. Quasar biasing 

The quasar biasing could be fairly complex in principle, but again in the linear regime that we consider here, one 
can safely assume the scale-independent and linear bias (e.g., Matsubara 1999): 

£ QS o(r,*) = [b QSO (z)] 2 U(r,z) = [b QSO (z)] 2 D(z) 2 U(r,0). (7) 

Then it should be emphasized here that all the above three effects (the light cone effect, the rcdshift-space distortion, 
and the quasar biasing) arc linear in a sense that equations (5) to (7) yield a linear relation between the quasar 2PCF 
in redshift space on the light-cone and the matter 2PCF in real space at z = 0. In particular the three effects simply 
change the amplitude of the matter 2PCF but none changes the positions of the baryon bump and the zero-crossing 
scale. This is valid only in the linear regime, and significantly simplifies the analysis. Consequently we can simply 
introduce an effective biasing factor, b e ff. 

£ Q S0(s)=&cff 2 £mM), (8) 

and 



b, 2 



eff = 



dz 



dVA- 1 ( dN QSO (z) 



dz J \ dz 



We note a possible systematic effect that gives rise to scale-dependence in the effective bias: if mean redshifts of 
quasar pairs which contribute to small-scale and large-scale correlations are significantly different, the biasing measured 
on the light-cone may acquire artificial scale-dependence. In other words, the redshift dependence of the selection 
function combined with that of the scale-independent bias may induce additional scale-dependence in the effective bias. 
For instance, Tegmark et al. (2004) discussed in detail how the luminosity dependence yields such scale-dependence of 
biasing using a subset of the DR2 SDSS galaxy catalog. In the current sample, however, this is not a serious problem 
for two reasons. First, while the luminosity and color dependence of galaxy clustering is well established (Kayo et al. 
2004; Zchavi et al. 2004), it is not yet clear that quasars also exhibit similar strong dependence on their intrinsic 
properties. Croom et al. (2004), for instance, find no significant luminosity dependence in the clustering of 2dF QSOs. 
Second, as is clear from the QSO distribution plot and dN/dz (Figs. 1 and 2), all separation bins of the 2PCF that 
we use in the analysis are not dominated by those in any specific redshift range, rather they are contributed equally 
over the entire redshift range. To be more specific, we compute the mean redshifts of pairs, z, at separation bins of s 
assuming fi m = 0.3 and ft A = 0.7. We find (z, s) = (1.40, lS/^Mpe), (1.41, 25/i" 1 Mpc), (1.41, lOO/i^Mpc), (1.41, 
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200/i" 1 Mpc), (1.41, 400/i~ 1 Mpc), (1.42, 800/i _1 Mpc), and (1.41,1600/i _1 Mpc). This reasonably justifies the fact that 
we neglect the luminosity dependence of quasar biasing. 

Thus wc proceed as follows; first we specify the values of fib-, A , and h, then compute £ m (r, 0). Since we assume 
the spatially flat model, O m is given as 1 — A . As discussed below, wc do not have to specify the value of as, the 
top-hat mass fluctuation amplitude at 8/i _1 Mpc, as long as we do not intend to solve for &qso(z). Second we evaluate 
£qso( s ) from the quasar sample using the same values of O m , A , and h. Finally, we compute the \ 2 of the observed 
and the model 2PCFs by varying the value of & e ff- In doing so, we use both diagonal and off-diagonal elements of the 
covariance matrix: 

X 2 = X] (^QSO,obs,i - &eff£m,i) (^QSCXobsj ~ &eff£m,j) Cy " 1 , (10) 



where C%f~ is the inverse matrix of the covariance matrix and the subscripts, i and j, are indices of separation bins. 
We repeat the above procedure for < Ob < 1 an d < A < 1, and derive constraints on the two parameters. Wc 
adopt h = 0.7 as a fiducial value, but also show results for h = 0.6 and h = 0.8 as well, although the /i-dependence is 
very weak. 

The resulting constraints are independent of the value of as and also of the evolution model of 6qso( z )- Nevertheless 
it is also interesting to estimate the quasar biasing factor even in a model-dependent manner. We adopt a simple 
evolution model by Fry (1996): 



This model neglects the possible luminosity /color dependence of quasar biasing and the finite lifetime of quasars, and 
may need to be improved to some extent; we note that Croom et al. (2004) find no significant luminosity dependence 
in the clustering of 2dF QSOs. Nevertheless equation (11) gives us a useful measure of the expected degree of the 
quasar spatial biasing at z = 0. 

4. Results and Discussion 

In performing the statistical analysis, we recompute £Qso,obs(s) for different values of Oa (0 < A < 1 at an interval 
of AO A = 0.05). Figure 8 shows the observed 2PCFs of the SDSS quasars for several different values of A (the n= 1 
binning offset). Each data point and associated errorbar is slightly shifted horizontally for clarity. Unlike in Figures 
5 to 7, the baryon bump is not directly visible due to the large bin size. Because of the statistical limitation, we also 
do not directly detect the anti-correlation that is predicted beyond a scale of 200/i~ 1 Mpc (in the case of fib = 0.04 
and A = 0.7 universe). However, we find that there is a break in the slope of the 2PCF around 100/i _1 Mpc. While 
?QSO,obs(s) is always positive on the scales less than 100 /i~ 1 Mpc, it is not inconsistent to be negative on the larger 
scale for any sets of cosmological parameters that we surveyed here. Thus this break may be interpreted to put a 
lower limit on the zero-crossing scale of CQSO.obs(s)- 

Figure 9 shows the 2PCFs of SDSS quasars using the 6 different binning offsets (O m = 0.3 and A = 0.7). The break 
of ^Qso,obs(s) around 100/i _1 Mpc is quite robust against the different choice of the binning. Since the data at different 
bins are strongly correlated, we do not use this entire dataset for our analysis; this is just to show that the binning 
makes little difference. In what follows, we perform statistical analysis for the n = 1 and n = 4 offsets, separately. 

Figure 10 shows the corresponding contours of 



where fo c ff , m in is the value of 6 e ff which minimizes x 2 f° r a given set of fib and A , and Xmin ^ s the global minimum 
value in fib, flA and b c g parameter space. The three curves represent the 68%, 95% and 99.7% confidence levels. 
Figure 11 is the same plot but employs the n = 4 binning offset (see Fig. 9). 

The cosmological parameters inferred from the WMAP results (Spergel et al. 2003) are in agreement with our 
constraints in Figures 10 and 11 at the 68% confidence level. Given these uncertainties, therefore, we express our 
constraints as allowed regions of fib and fl\. For this purpose, we fit the mean and 2a contours as straight lines and 
obtain 




bQSo(z) = 1 + 



&qso(z = 0)-1 
D(z) 



(11) 



Ax 2 (flb,flA) = X 2 (O b ,0 A ,6 offjm i n ) -Xmi 



(12) 



2.10b + A = 0.89 (mean), 
0.80 - 2.80b < A < 0.90- 1.40b (2a) 
for the n = 1 binning offset (Fig. 10), and 
3.60b + A = 0.99 (mean), 
0.69 - 5.10b < A < 1.0 -2.30b (2a) 



(13) 
(14) 



(15) 
(16) 



No. 
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for the n = 4 binning offset (Fig. 11). These fitted lines are also plotted in Figures 10 and 11. The differences of 
the slopes of contours between the n = 1 and n = 4 binning offsets mainly come from the smallest-scale behavior of 
^Qso.obs(s), but the resulting constraints arc fairly insensitive to the different binning offsets (compare Figs. 10 and 
11). The dependence of our constraints on the Hubble constant is illustrated in Figure 12. The upper and lower 
panels are derived by employing the n = 1 and n = 4 offsets, respectively. The left and right panels show the results for 
h = 0.6, and h = 0.8, and the central panels plot the results for h = 0.7 (identical to Figures 10 and 11) for comparison. 
Except for the small difference, the overall conclusion is that our constraints are very insensitive to the choice of h. 

Finally, let us consider the corresponding biasing parameter of quasars implied from the current study. Unfortunately 
the current sample cannot constrain the redshift-dependence on the bias evolution. Thus we have to specify the value 
of as (at z=0) and the evolution model of quasar biasing as well, although all the constraints presented above are 
independent of those assumptions . For simplicity and clarity, we adopt = 0.04 and Q\ = 0.70. Then we find that 
6 G ff is 1.13<78 -1 f° r 0.6 < a < 1.0. If we define a "mean" quasar bias as 



dVc 
dz 



dN QSO (z) 
dz 



l + |?(z) + i/3 2 (z) 



dz 



dV c 

dz 



1 f dN QSO (z) 
\ dz 



D 2 (z) 



(17) 



we find that b mean = 1.59<78 -1 (the mean redshift of the sample is fisl.46). This is a measure of the degree of quasar 
clustering if the quasar biasing is independent of redshift. We think this result is consistent with Croom et al. 2004. 
If we adopt another bias evolution (eq.[ll]) instead, the current quasar sample is fitted to have 6Qso,Fry(z = 0) = 
0.54 + 0.83CT8- 1 . 



5. Summary 



We have presented the first result of the SDSS quasar clustering analysis. In particular we have focused on the 
behavior of the quasar two-point correlation functions (2PCF) at large scales, and explored the embedded baryonic 
signature through the overall shape of the quasar 2PCF. The measurement is still statistically limited, but we found 
that the zero-crossing scale of the quasar 2PCFs is larger than 100/i -1 Mpc. The estimated quasar 2PCFs are in good 
agreement with those predicted in linear theory combining the light-cone effect, the rcdshift-space distortion and the 
scale-independent linear biasing model. 

To proceed further, we have attempted to derive joint constraints on the baryon density parameter f2b and the 
cosmological constant J7a under the assumption of the spatially flat universe. Our constraints are approximately 
expressed by the relation: 0.80 — 2.8f2b < < 0.90 — 1.4f2b at the 2<r confidence level (n = 1 binning offset). While 
the statistical significance is not competitive, the constraint is consistent with the independent results from other 
observations such as WMAP (Spergel et al. 2003) within the la level. 

The current analysis can be improved in many different ways; to consider non-flat + Oa ^ 1) models, to explore 
a more general equation of state of the universe than the cosmological constant (e.g., Yamamoto 2004), to consider 
the anisotropy of the 2-dimcnsional correlation function (Matsubara, Suto 1996), to incorporate the gravitational 
nonlinearity effect which is not important on ~ 10/i -1 Mpc scales but still cannot be entirely neglected (Changbom 
Park, private communication), to take into account the finite angular size of the pair-separation (i.e., without adopting 
the distant-observer approximation; see Matsubara 2000), and to elaborate on the statistical modeling, in particular, 
the error estimation. We plan to perform such improved analyses in due course. 

The wedge diagram of the current quasar sample, Figure 2, clearly illustrates the isotropy and homogeneity of the 
universe viewed over large scales, in marked contrast to the complex clustering pattern of the local universe traced by 
galaxies (e.g., Fig. 3 in Hikage et al. 2003). Therefore the most impressive lesson that we have learned from the current 
analysis may be the fact that we can indeed decipher cosmological information on f2b and Oa which is imprinted in 
the apparently almost homogeneous distribution of quasars at such high redshifts. 
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Large Scale Clustering of SDSS QSOs 




Fig. 1. Equatorial coordinate distribution of SDSS quasars that we use in the current analysis (19986 in total). 




M^mag] 



2. Wedge diagram of the quasar distribution (19986 quasars in total for 0.72 < z < 2.24) corresponding to Figure 1. 
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Fig. 3. Scatter plot of redshift and i-band absolute magnitude (Upper panel) and the corresponding redshift his- 
togram (Lower panel) of our quasar sample (0.16 < z < 2.24 and 15.0 < mi^ c < 19.1). The width of the red- 
shift bin is 0.08. To compute the absolute magnitude, we assume that the energy spectrum of quasars fol- 
lows a single power-law with a spectral index of —0.5, and adopt the cosmology of f2 m = 0.3 and = 0-7. 
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Fig. 4. The two-point correlation functions of our samples, BEST DR3 and TARGET DR3 (the n = 1 binning offset) for 
Qm = 0.3, Qa = 0.7 and h = 0.7. Filled circles and triangles correspond to our samples (DR4 and DR3), and filled squares 
and stars correspond to the BEST and TARGET DR3 quasar samples in the catalog by Schneider et al. (2005). 
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Fig. 5. Matter two-point correlation functions in linear theory (in real space) at z = for different values of fl^,. The CDM 
transfer function of Eisenstein, Hu (1998) is adopted. We assume the spatially flat model (f2b + fid + £1a = 1), and h = 0.7. 




Fig. 6. Same as Figure 5 but for different values of Qa with H^, = 0.04 and h = 0.7. 




Fig. 7. Same as Figures 5 and 6 but for different values of h with = 0.04, = 0.23, and Qa = 0. 
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Fig. 8. The two-point correlation functions of SDSS quasars (the n = 1 binning offset) for several different values of Qa (h 
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Fig. 9. The binning offset dependence of the two-point correlation functions of quasars. Results for the 

six different binning offsets (see eq.[2]) are plotted. We adopt Q m = 0.3, Qjy = 0.7, and h = 0.7. 
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Fig. 10. Contour plot of A\ 2 (eq.[12]) on the - Qa plane (h = 0.7). The three contour curves represent the 68%, 95% 
and 99.7% confidence levels. The cross and the star indicate the values preferred by WMAP (Ob = 0.04 and Q\ = 0.73), 
and our best-fit model (Qt> = 0.1 and Qa = 0.7). We adopt the n = 1 binning offset in equation (2). The solid and 
dotted straight lines indicate the linear fits to the mean and the 2<r contour levels (eqs.[13] and [14]), respectively. 




Fig. 11. Same as Figure 10 but for the n = 4 binning offset. The soiid and dotted straight lines 

indicate the linear fits to the mean and the 2tr contour levels (eqs.[15] and [16]), respectively. 




Fig. 12. Same as Figures 10 and 11 but for h = 0.6 {Left panels), h = 0.7 (Center panels), and 
h = 0.8 (Right panels). Upper panels: n = 1 and Lower panels: n = 4 binning offsets. 



